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ABSTRACT 



X 



Context. One of the most intriguing features discovered by Swift is a plateau phase in the X-ray flux decay of about 70% of the 

afterglows of gamma-ray bursts (GRBs). The physical origin of this feature is still being debated. 

Aims. We constrain the proposed interpretations, based on the intrinsic temporal properties of the plateau phase. 

Methods. We selected and analyzed all the SwiftfXKt GRB afterglows at known redshift observed between March 2005 and June 

2008 featuring a shallow decay phase in their X-ray lightcurves. 

Results. For our sample of 21 GRBs we find an anticorrelation of the logarithm of the duration of the shallow phase with redshift, 
with a Spearman rank-order correlation coefficient of r = -0.4 and a null hypothesis probability of 5%. When we correct the durations 
for cosmological dilation, the anticorrelation strenghtens, with r = -0.6 and a null hypothesis probability of 0.4%. Considering only 
those GRBs in our sample that have a well-measured burst peak energy (8 out of 21), we find an anticorrelation between the energy 
of the burst and the shallow phase duration, with r=-0.80 and a null hypothesis probability of 1.8%. 

Conclusions. If the burst energy anticorrelation with the shallow phase duration is real, then the dependence of the shallow phase 
on redshift could be the result of a selection effect, since on average high-redshift bursts with lower energies and longer plateaus 
would be missed. A burst energy anticorrelation with the shallow phase duration would be expected if the end of the plateau arises 
from a collimated outflow. Alternative scenarios are briefly discussed involving a possible cosmological evolution of the mechanism 
responsible for the X-ray shallow decay. 
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1. Introduction 

In the pie-Swift era the X-ray afterglows of gamma-ray bursts 
could be observed only after many hours after the burst, when the 
flux typically showed a smooth power law-like decay, f with 
an index of about a ~ 1. Hereafter, we refer to this as the stan- 
dard X-ray afterglow decay phase. The Swift mission (Gehrels 
et al. 2004) has revolutionized GRB studies in many respects by 
observing the X-ray afterglow phase from a few dozen seconds 
after the burst (e.g. Zhang et al. 2007). 

The shallow decay phase observed in the X-ray flux of about 
~ 70% of the afterglows (e.g. Panaitescu 2007) is one of the 
most intriguing features discovered by Swift. This phase usually 
becomes visible a few hundred seconds after the burst, after the 
steep decay in the prompt emission, and it lasts for ~ 1 - 10 
ks (Nousek et al. 2006; Zhang et al. 2006). No spectral evolu- 
tion is observed in either the 0.3-10 keV range during the shal- 
low phase or in the subsequent decay phase (e.g. Liang et al. 
2007; Butler and Kochevski 2007). This lack of X-ray spec- 
tral variations has suggested that the observed X-ray temporal 
steepening is not associated with the crossing of a characteristic 
synchrotron frequency (e.g. the cooling frequency). Optical af- 
terglow lightcurves often show a different behavior from those 
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in the X-rays (Liang et al. 2007). This idicates either that the 
X-ray and optical afterglow have different origins or that the mi- 
crophysical parameters determining the instantaneous energy in 
the electrons and magnetic field evolve in time (Panaitescu et 
al. 2006). Some of the models proposed to interpret the physical 
origin of the shallow decay phase are summarized in Zhang et 
al. (2007). 

Several studies have addressed the intrinsic properties of the 
X-ray shallow phase, in particular by testing whether a depen- 
dence exists between the intrinsic duration of the shallow phase 
and the burst energetics. Results obtained so far are discordant. 
An anticorrelation between the intrinsic duration of the shallow 
phase and the burst energetics has been found in some works 
(e.g. Sato et al. 2007; Dado et al. 2008), while in some others 
it was not (e.g. Liang et al. 2007; Nava et al. 2007). Sato et al. 
(2007) argue that these discrepancies may be associated with the 
large uncertainties affecting the burst energetics estimates and/or 
in modeling the X-ray shallow phase and estimating the tempo- 
ral break between the shallow and the standard phases. Some of 
the discrepancies might also be ascribed to the different size and 
quality of the GRB samples used by different authors. 

A linear dependence between the logarithm of the duration 
of the shallow phase, and the logarithm of the burst istropic 
equivalent energy would be expected if the temporal break be- 
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Fig. 1. X-ray afterglow lightcurves from Evans et al. (2007) 
where the temporal axis has been rescaled to the burst rest frame 
for 3 GRBs of our 'golden sample' a different redshift (see §2). 



tween the shallow phase and the standard phase were interpreted 
as 'jet break' time (tj) and the GRB energy corrected for beam- 
ing factor were constant (e.g. Frail et al. 2001). Indeed, the jet 
opening angle can be estimated from the time tj at which the rel- 
ativistic beaming (l/r(tj) where r(tj) is the fireball Lorentz fac- 
tor) becomes equal to the geometric beaming of the fireball of 
half-opening angle Oj, that is, 0j oc (Pj/Eisoy/^ (Sari et al. 1999), 
where £■,,„ is the equivalent isotropic energy. At that time (tj), the 
afterglow lightcurve decay steepens. The beamed corrected en- 
ergy is Ey - Eiso(l -cosOj) ~ EjsoOj. If Ey is constant, it follows 

that Eiso °^ fj'. The correlation found by Ghirlanda et al. (2004) 
between the intrinsic peak energies and the beaming-corrected 
burst energies tells that the relation Eiso °^ is still (nearly) 
valid for GRBs with similar intrinsic peak energy. In the jet sce- 
nario, the lightcurve steepening is expected to be achromatic. 
This condition is barely satisfied if we consider both the X-ray 
and the optical energy domains, since, as mentioned above, sev- 
eral X-ray shallow phases are not tracked in the optical regime. 
However, by restricting the energy range to the X-rays, the con- 
dition is satisfied since the lack of any spectral variation is a 
characteristic feature of the X-ray shallow phase and the subse- 
quent standard decay phase. 

In the present paper we consider a sample of GRBs with a 
well-monitored X-ray lightcurve and known redshift, which un- 
ambigously showed a shallow decay phase. We find clear evi- 
dence of a redshift dependence of the duration of the shallow 
decay phase. 



2. The sample and data analysis 

The X-ray afterglow lightcurves supplied by the UK Swift 
Science Data Centre at the University of Leicester were used 
(Evans et al. 2007). The sample was built by selecting in the 
Swift archive all GRBs observed in the period from March 2005 
to June 2008 with the following characteristics: i) a 0.3-10 keV 
XRT lightcurve featuring shallow behavior with temporal index 
ai < 0.8 (i.e. shallower than the 'standard' fireball model pre- 
dictions, e.g. Sari et al. 1998), over a temporal interval greater 
than 0.5 ks (so that the power-law decay index can be mea- 
sured accurately); ii) the shallow decay interval should not be 




dominated by features such as spikes or flares that may affect 
the measurement of the decay index; iii) since our analysis con- 
centrates on the study of the intrinsic properties of the shallow 
phase, we considered only GRBs at known distances. In this way 
we selected 21 GRBs (out of a total of about 60 bursts at known 
redshift for which the statistics were good enough to carry out a 
detailed temporal analysis). The sample, as well the redshift of 
each burst, is given in Table 1 . In Figure[T]the lightcurves of three 
GRBs from our 'golden sample' (see below), taken from Evans 
et al. (2007), are plotted after rescaling of the temporal axis to 
their rest frame. The shape of the lightcurves clearly shows the 
well known 'canonical' behavior (e.g. Nousek et al. 2005). Since 
our goal is to measure the duration of the shallow phase, we fit- 
ted the observed lightcurves with power laws as: 



(1) 



In Eq. (1), f < fi indicates the interval where the typical 
initial steep decay from the prompt is observed, t[ < t < t2 
corresponds to the shallow decay phase, and f > f2 to the subse- 
quent standard decay. In this work, we consider t2 as a measure 
of the duration of the shallow decay. For GRB 060614 and GRB 
060814, we excluded the initial steep decay from the fit because 
it could not be well-fitted by a power-law decay. We also ex- 
cluded from the fit those intervals in which flares were present 
for f < fi in the lightcurves of other GRBs (see last column of 
Tab. 1). 

Table [T] shows the selected dataset where f2 represents the 
observed epoch at which the shallow phase steepens to the stan- 
dard decay (calculated from the burst onset as determined by 
Swift/BAT), and ai represents the temporal index of the shal- 
low decay region. Errors are given at the Icr confidence level. 
In some cases Eq. (1) provided a poor approximation of the 
steepening from the shallow phase to the standard one. We thus 
checked whether other estimates of the temporal break between 
the shallow phase and the standard phase obtained assuming 
more complex models (e.g. Willingale et al. 2007; Ghisellini et 
al. 2008) provided different results and we find no significant 
differences within the uncertainties, except for two case^. 

We first checked whether any common intrinsic value t'^ of 
the epoch at which the X-ray shallow decay ends (t^ = f2/(l +z)) 
exists for all GRBs. The observed epoch f2 covers 3 orders of 
magnitudes (0.5-80 ks). The intrinsic epoch still covers a wide 
range of values (0.2-50 ks. Tab. 1). We find no evidence of clus- 
tering around any particular value. We then checked whether 
there is any redshift dependence on f2. We found that an anti- 
correlation exists between the logarithm of t2 and z, with a rank 
correlation factor of -0.4. With 19 degrees of freedom, the null 
hypothesis is rejected at a 95% confidence level. After correcting 
t2 for cosmological dilation, we found that the anticorrelation is 
strengthened, with a rank correlation factor of -0.6: the null hy- 
pothesis is now rejected at 99.6% (non-directional probability) 
confidence level (Fig.|2]l. That the anticorrelation already present 
for t2 becomes more significant after correction for cosmological 
dilation provides evidence that the correlation is genuine and not 
biased by the redshift correction. To confirm this result, we se- 
lected a "golden sample" of 10 GRBs from our original 21 GRB 
sample by considering only those GRBs with the best XRT cov- 
erage in all the three typical SwiftpCRT X-ray lightcurve compo- 



' Willingale et al. (2007) estimated for GRB 060607A a duration of 
r„ = 56*4. ks and Ghisellini et al. (2008) for GRB 050319 estimated 
= 7 ks (the latter is in the rest frame). 
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Table 1. Selected sample of GRBs with known redshift that presents a shallow decay phase. 

From Amati et al. 2008; From Nava et al. 2007. 



nents (e.g. Nousek et al. 2005), which are an initial steep decay 
followed by the shallow and then standard decay. Despite the 
decrease in the sample of GRBs, the anticorrelation between the 
logarithm of t2 and z persists, with a rank correlation factor of - 
0.85. (The null hypothesis is rejected at 99.2% confidence level.) 
The golden sample is marked with red open circles in Figure |2] 
and in boldface in Table 1 . 

These findings cannot be interpreted as due to an energy de- 
pendence of the duration of the shallow phase as one may con- 
clude for example from Figure [T] where the shallow phase at 
energies of 1-37 keV (for ORB 060714 at z=2.711) in the rest 
frame is shorter than the one at 0.5-15 keV (for GRB 060729 
at z=0.54). In fact, it is well known that the hardness ratio, de- 
fined as the flux ratio in the 0.3-1.5 keV and 1.5-10 keV energy 
bands, does not show any evidence of variations along the shal- 
low phase (e.g. Butler and Kocevski 2007; Liang et al. 2007). 

The anticorrelation of the intrinsic duration of the X-ray 
shallow phase with redshift that we discussed above may be 
a consequence of an anticorrelation of f, with burst energy. 
Indeed, the GRBs that we observe at high-redshift are on av- 
erage more energetic than the low-redshift ones due to a simple 
selection effect. We cannot see faint GRBs at large distances. 
Therefore, shorter plateaus might be observed more frequently 
at high redshift because associated with more energetic GRBs. 
Alternatively, the anticorrelation of t'^ with z can be explained as 
a cosmological evolution of the mechanism that gives rise to the 
shallow decay. We briefly discuss both possibilities here. 

2.1. The burst energy dependence of ^ 

To verify the dependence of the energy of the burst from is 
not an easy goal since the GRB energetics, usually estimated by 
Eiso, is often affected by large uncertainties in the burst spec- 
tral parameters, the peak energy Epeak of the EFe spectrum in 
particular. For this reason, we considered only those GRBs of 
our sample for which precise measurements of Epeak are avail- 
able (e.g. Amati et al. 2008; Nava et al. 2007). For 8 GRBs of 
our sample (3 of which are part of the 'golden sample', see Tab. 



1 and Figure |2]l, we find that t'^ anticorrelates with Eiso, with a 
rank correlation factor of -0.80 and null hypothesis rejected at 
98.2% confidence level (Fig.|2]i. If we restrict ourselves to con- 
sidering only those GRBs with similar intrinsic peak energy, we 
still find marginal evidence of the anticorrelation, although the 
statistics are poor (Fig.|2l). Even though a firm conclusion could 
not be reached, this result is consistent with the results obtained 
by Dado et al. (2008) and Sato et al. (2007), while at odds with 
findings by Liang et al. (2007) and Nava et al. (2007). 

As mentioned in § 1 , the anticorrelation between and Eiso 
suggests that t2 could be considered as the jet break time of (the 
part of) the fireball that gives rise to the X-ray afterglow (see § 1). 
Several suggestions have been made to reconcile the jet interpre- 
tation with features observed in some GRB, such as a chromatic 
evolution of the break or the temporal break at later epochs. For 
example, the two-component jet model (e.g. Peng et al. 2005; 
Racusin et al. 2007) explains the lack of a simultaneous optical 
break as a deficit in the optical emission from the narrower of 
the two jets, responsible for the two X-ray temporal breaks. In 
another model, Ghisellini et al. (2007) propose to interpret the 
observed X-ray plateau as the sum of two components: the late 
prompt emission (internal shock between late emitted shells) and 
the afterglow. The temporal break at the end of the X-ray shallow 
decay phase is the proof that the Lorentz factor of the late shells 
(typically smaller than that of the external shell giving rise to 
the afterglow) has reached the l/9j value (see §1). Whether the 
break is tracked in the optical band depends on the relative in- 
tensity of each component. The second break at later times is 
produced when the Lorentz factor of the shell producing the af- 
terglow has reached the l/Oj value. 

2.2. Evidence of cosmic evolution? 

Given the uncertainties affecting the burst energetic evaluations 
and the poor statistics available to definitively confirm the anti- 
correlation between Eiso and 4' speculate on another possible 
interpretation of the redshift anticorrelation with where the 
shallow phase depends on an external component that evolves 
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3. Conclusions 
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In this work we analyzed 21 GRBs with known redshift that 
feature a shallow phase in the X-ray lightcurve. Our main re- 
sult is a clear anticorrelation of the intrinsic duration of the X- 
ray shallow phase with redshift. Considering only those GRBs 
in our sample that have well-measured burst peak energy, we 
find marginal evidence for burst energy anticorrelation with the 
shallow phase duration. The latter anticorrelation would be ex- 
pected if the observed temporal break (fi) arises from a col- 
limated outflow. In this case, the fj anticorrelation with z can 
be interpreted as the evidence of a selection effect since high- 
redshift bursts with lower energies and shorter plateaus would 
be missed. However, a larger sample of bursts at known red- 
shift with well-measured burst spectral parameters is required to 
definitively assess the Ejso anticorrelation with t^. In an alterna- 
tive scenario, the shallow phase may arise from a mechanism 
that operates differently at high redshift, such as for example 
from X-ray dust scattering or an evolution of the intrinsic GRB 
properties as for example the inner engine (Guetta et al. in prepa- 
ration). Finally we note that by virtue of its redshift dependence 
(though with large scatter), the observed X-ray shallow phase 
duration {to) may be regarded as an additional figure of merit 
to single out high-redshift GRBs directly from X-rays observa- 
tions. This might provide useful information for burst follow-up 
campaigns at optical and NIR wavelengths. 
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Fig. 2. Top panel: redshift versus the intrinsic duration of the 
shallow phase. The observed anticorrelation has a null hypothe- 
sis rejected at > 99.6% confidence level. The filled black circles 
are the GRBs with measured peak energy and those encircled in 
red represent the 'golden sample' (in bold face in Tab. 1). The 
rest of the sample is plotted with open circles. Bottom panel: the 
isotropic equivalent energy of the burst versus the intrinsic du- 
ration of the shallow phase. The observed anticorrelation has a 
null hypothesis rejected at > 98.2% confidence level. Black cir- 
cles, red triangles, and blue square indicate GRBs with similar 
intrinsic peak energy (400 < Epgakj ^ 600 keV, Ep^akj < 200 
keV and Epeakj > 1000 keV respectively, see Table 1. 



with redshift. This is the case, for example, if the shallow phase 
is produced by the interaction of the X-ray emission with the sur- 
rounding interstellar dust (e.g. Klose 1998,1999; Shao and Dai 
2007,2008). The expected decrease in the interstellar dust con- 
tent with redshift might be the reason for the observed anticor- 
relation between the duration of the shallow phase and redshift. 
Moreover, it is expected within this model that the X-ray and op- 
tical temporal breaks are in general uncorrelated, in agreement 
with a number well-sampled afterglow lightcurves (Liang et al. 
2007). However, this interpretation faces problems in explaining 
the lack of the predicted spectral variation in the X-ray spectra, 
as akeady pointed out by Shen et al. (2008). A possible solu- 
tion might involve complex dust distribution along the line of 
sight. Other scenarios are still possible, as for example if the X- 
ray shallow phase depends on the intrinsic GRBs properties (e.g. 
inner engine) that may evolve within cosmological time scales. 
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